UMD-DOE/40762-516 
DAMTP-2012-20 



Neutrino Emission from Helium White Dwarfs with Condensed Cores 

Paulo F. Bedaque* and Evan Berkowitz^ 

Maryland Center for Fundamental Physics, 

Department of Physics, 

University of Maryland, College Park, MD USA 



(N 

o 



Aleksey Cherman-'- 

Department of Applied Mathematics and Theoretical Physics, 

University of Cambridge, Cambridge CB3 OWA, UK 

The possibility that ions in a helium white dwarf star are in a Bose-Einstein condensed state has 
been explored recently. In particular, it has been argued that the resulting novel quantum liquid 
has a new kind of quasiparticle excitation with a phonon-like dispersion relation. We investigate 
the neutrino emission rate due to this gapless state and the resulting impact on the total luminosity 
of helium white dwarf stars, as a possible observable way of detecting this exotic phase. If the 
condensation temperature for the quantum liquid state, which is currently not known very precisely, 
turns out to be high enough, our calculations indicate that neutrino emission due to the gapless 
mode would make a large contribution to the total luminosity of the helium white dwarf stars. 
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Introduction. Stars which are not too heavy become 
white dwarfs upon running out of fuel for nuclear fu- 
sion. White dwarf stars evolve by cooling, and have two 
cooling mechanisms. Very young white dwarfs lose en- 
ergy mostly through the decay of massive plasmon quasi- 
particles to neutrinos for temperatures ^ 10^ Kelvin [1- 
3] . At these temperatures the effect of neutrino emission 
may be detected in frequency changes of non-radial os- 
cillations of the star [4, 5]. However, plasmons become 
Boltzmann-suppressed for T < 10^ Kelvin, and below 
this temperature the dominant energy loss mechanism 
for white dwarfs is believed to be electromagnetic radia- 
tion from the surface of the star. 

Here we discuss the neutrino emission rate for helium- 
core white dwarfs (He WDs). He WDs are much rarer 
than the usual carbon-oxygen white dwarfs (CO WDs), 
and are only formed if a red-giant stage star in e.g. a bi- 
nary system loses much of its hydrogen envelope before 
helium burning can begin. The first He WDs candidates 
were detected only recently [6-8] in the globular cluster 
NGC 6397. In Ref. [8] it was observed that the He WD 
sequence terminates early compared to the CO WD se- 
quence, so that there is an apparent dearth of faint He 
WDs. This unexpected result motivates a reexamination 
of the physics of super-dense helium plasmas. 

Recently it was realized that for low enough tempera- 
tures, the helium nuclei in a high density helium plasma 
should Bose-condense, forming a 'nuclear condensate' [9- 
12]. At the relevant temperatures and densities, the 
specific heat c„ of the condensed phase turns out to 
be ~ 10~^ times smaller than the c^ of the classical 
plasma phase, and ^ 10^^ times smaller than the spe- 
cific heat of an ion lattice, if one were formed. This 
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implies that He WDs with condensed core should cool by 
electromagnetic radiation significantly faster than WDs 
with uncondensed cores. Whether this cooling speedup 
is large enough to be observable depends on detailed 
modeling[10, 13], and the situation is currently not en- 
tirely clear. 

In Ref. [12] we showed that the condensed phase of 
super-dense helium has a (naively) unexpected gapless 
quasiparticle mode, which has important implications for 
the physical properties of nuclear condensates. This low- 
temperature mode is a collective excitation of the elec- 
trons and (condensed) ions. The mode is a density os- 
cillation with a linear dispersion relation, and so it is a 
sound-like mode which propagates at zero temperature. 
However, it is not the zero sound familiar from Fermi 
liquids, and for lack of a better name we will call this 
gapless mode 'half-sound'. In Ref. [12] we examined the 
effect of the half-sound on the specific heat of the nu- 
clear condensate, which helps set the rate of cooling via 
electromagnetic radiation. However, the existence of the 
gapless half-sound mode should also affect the neutrino 
cooling rate, since it is not Boltzmann suppressed for any 
temperature T, in contrast to the gapped plasmon mode. 
The purpose of this paper is to investigate whether the 
annihilation of half-sound quasiparticles into neutrinos 
might make an appreciable contribution to the total lu- 
minosity of He WDs, and hence to their cooling rate. 

Ignoring quantum effects, at zero temperature and at 
the relevant densities the helium nuclei would arrange 
themselves, due to the Coulomb repulsion, into a lat- 
tice immersed in a neutralizing sea of electrons. At a 
higher temperature Tmcit this lattice melts. We can es- 
timate Tmcit by equating the typical Coulomb energy to 
the temperature; the resulting estimate is Tmoit ~ an^'^ , 
where n is the density. On the other hand, since helium 
ions are bosons, at low temperatures they would be ex- 
pected to Bose-condense once their thermal de Broglie 
wavelength exceeds the inter-ion separation. The tem- 
perature at which such an ion condensate would melt 



can be estimated by equating the thermal wavelength 
of the ions to the inter-ion distance, leading to the es- 
timate Tcond ^ v?/^ /M. This shows that at sufficiently 
high density, there is a range of temperatures where the 
ions will be found in a condensate, and not in a clas- 
sical plasma or a lattice as one might have expected. 
For plasmas made from heavier elements this quantum 
liquid region would open up only at unrealistically high 
densities, but for a helium plasma the relevant densities 
are expected to be reachable in He WDs. The melt- 
ing temperature Tmeit has been analyzed carefully in the 
past and is well known[14-16], at least for the param- 
eters relevant for CO WDs. Tcond, however, is less well 
known. It is expected to be higher than the free gas result 
^cond ~ 1-49 T?!"^ jM on account of the strong repulsion 
between ions[17]. This expectation is born out by the 
computation in Ref. [18], which found a large enhance- 
ment over T'Jojj'Jj of about an order of magnitude; however 
the effects of the half-sound mode were not included in 
Ref. [18]. With the estimates of Ref. [11, 18], the region 
in parameter space where a condensate is expected to oc- 
cur overlaps with the relevant densities and temperatures 
in white dwarfs. 

Effective field theory. The coupled system of non- 
relativistic helium ions, electrons, and photons turns 
out to support three quasiparticle modes once the ions 
condense [12]. The two gapped modes are the plas- 
mon mode and the transverse photons, with gaps set 
by LUp — ATrZa/rrired £^nd jti^ — AnZ'^av'^/M respec- 
tively, where Z = 2 is the charge of the nuclei, M is 
their mass, a ~ 1/137 is the fine structure constant, and 
"rrired = t^eM / {Z fjLe + M) is the reduced mass of the ion- 
electron system, where /Xg is the electron chemical poten- 
tial. The chemical potential is related to the ion density 
n via /Zg = ml + kp^kp/Sn'^ = Zn, where rUe is the 
electron mass. At typical white dwarf conditions (cen- 
tral density of 5 x 10^^ g/cm^ and T < 3 x 10^ K ~ ttia 
), the gapped modes are Boltzmann suppressed, and the 
low-energy physics can be described in terms of the in- 
teractions of the gapless half sound mode, which has the 
Lagrangian 
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Ch = -H [-d^a + c'hV'] H 



(1) 



where ch — rnA/'nis, where rris — Aapiekp /t: is the De- 
bye screening mass due to the electrons, and the ellipsis 
denote interaction terms. Numerically, for the previously 
mentioned representative white dwarf, ch ^ 0.001c, 
where c is the speed of light, i.e., ch ~ 300km/s. 

Since iJ is is a gapless mode, the decay of H to neu- 
trinos is kinematically forbidden, and the leading con- 
tribution to the neutrino emissivity is from the annihi- 
lation reaction HH — >■ uD. The neutrino emissivity (en- 
ergy emitted per time per volume) Q can be estimated as 
Q ~ (G|./M^)T^^. Here the dependence on Fermi con- 
stant Gf is the usual one for weak-interaction processes. 
Meanwhile, the dependence on M is a consequence of the 
fact that the relevant interaction term must be quadratic 
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FIG. 1. The annihilation of two H particles into a neutrino 
pair. The dashed lines are the incoming H particles, and the 
solid lines are the outgoing neutrino/antineutrino pair. 



in H . The M dependence then follows by noting that the 
H mode is derivatively coupled, so that each power of H 
comes with a power of spatial momentum p, but since 
the helium nuclei are non-relativistic, p must enter in the 
combination p^ /2M . What such naive estimates cannot 
tell us, however, is the dependence of Q on cr, which 
turns out to give a parametric enhancement. To see that 
an enhancement is possible, note that the number den- 
sity of H quasiparticles is a function of their energy chP, 
so temperatures should enter as (T/ch), which would 
yield an enhancement of order cjj^^ ~ 10^'^. However, 
this is clearly an overestimate, since the phase space for 
the decay vanishes as c — ?> 0, and a detailed calculation is 
necessary to determine how ch enters Q. We will demon- 
strate that the correct dependence is Q oc cjj^ , yielding 
a nontrivial enhancement. 

Matching to the SM. To compute the neutrino emis- 
sivity (energy loss to neutrinos per time per volume), we 
must derive the coupling between the half-sound mode 
and neutrinos, which we do by a matching calculation be- 
tween the Standard Model (SM) down to the half-sound 
EFT in Eq. (1). The relevant coupling between H and 
neutrinos will be mediated by the neutral Z bosons, since 
the Z can decay to directly to neutrinos which can then 
escape the core, in contrast to the charged W^ bosons 
that would decay into a Pauli-blocked electron and a neu- 
trino. 

The 77 is a collective mode arising from the collective 
behavior of electrons and helium nuclei, and in general 
the H field will couple to the weak sector through both 
its constituents. In the SM, Z's coupling to the leptons 
is directly supplied by the Standard Model. The relevant 
piece of the Lagrangian is 
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67^(1 - 75)e -I- sin^ Owej'^e 



where 5 is a weak coupling constant, Ow is the Wein- 
berg angle, the fields v and e are neutrinos and electrons 
respectively. 



The SM couplings of the Z to quarks are described by we see from Eq.(6) that the H-Z couphng is^. 
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wyf^ (l — j5)u~ -sin^ Owu'j'^u (3) 



-d-f^il - 75)d + - sin^ Bwdj^'d 



We now need to deduce how a hehuni nucleus interacts 
with the Z, which we do by a series of matching calcu- 
lations. First, let us consider how the Z couples to the 
nucleon isospin doublet N. The Z must interact as 



C 



Z-had. 
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where gA ~ 1-26 is the nucleon axial charge, As = 
—0.16 ± 0.15 the strange axial charge, the matrix T^ — 
T^/2 the third component of weak isospin, Q the elec- 
tric charge in units of the fundamental charge, and the 
dots indicating higher dimensional terms which we will 
neglect. Notice that the vector coupling is not renormal- 
ized as it is a conserved current. 

We now take the nonrelativistic reduction of (4) so that 
we can more easily read off how a condensate of nuclei 
should interact with Z . It is simple to verify that 
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ZqN^T^N - {gA + As)ZiN''a'T^N 



sin^ OwZqN^QN 
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So, we see that the Z couples to the total weak isospin. 
The helium nucleus has zero weak isospin, and thus does 
not couple to the Z through the first term at all. Be- 
cause the ground state of the ^He nucleus is dominated 
by configurations where the total spin of the four nucle- 
ons vanishes (only 10% — 15% of the wave function has 
non- vanishing spin [19]), we expect the second term to 
be negligible, and so we will ignore the axial coupling. 
We are left with 
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where the field V' represents the ^He nuclei. Note that 
this coupling is suppressed by the Weinberg angle. In 
a fictitious world where sm9w = 0, the Z would be 
the gauge boson related to the third component of weak 
isospin and thus would not couple at all to an isosinglet 
source like ^He. 

Finally, recalling from Ref. [12] that H is related to tp 
via 
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Because the regime we are interested in includes temper- 
atures up to roughly ^ lO^K « 10~^GeV, which is vastly 
smaller than the mass of the Z, mz ~ 90GeV, we can 
immediately integrate out the Z, directly connecting the 
H excitations with neutrinos. A simple calculation yields 
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where in the second line we used {g/ cos 9 
8Gi.A/|/\/2. 

Neutrino emissivity. We are now in a position to com- 
pute the neutrino emissivity Q, which can be written as 
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n(fco)n(fc^)(fco + fco)r(A:,fc')- (10) 



where n^ko) is the Bose-Einstein distribution (exp(/3fco) — 
1)~^, /3 = 1/T, fcg, /cq are the energies of the annihilating 
H modes with momenta of magnitude k, fc', and r(fc, fc') 
is the annihilation probability per time T per volume 
V. Writing the annihilation amplitude for the process 
of Fig. 1, which can be read from (9), the annihilation 
probability per time per volume evaluates to 
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and we are summing over the helicities and momenta p, p' 
of the neutrinos with energies po « p, pg « p' in the final 
state. The momentum integral can be evaluated in closed 
form, yielding 
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Plugging this result into (10), the remaining integral over 
the momenta of the annihilating H quasiparticles can be 
evaluated analytically in the approximation that c// <C 1, 
with the result 
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^ H mixes with Z but this weak force effect is neghgible. 
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FIG. 2. The "L^" (blue) lines are the electromagnetic lu- 
minosities from (15) and from the model in Ref. [20] (the 
dotted line is our extrapolation). The (yellow) "L^" line is 
the neutrino luminosity for a fictitious constant-density star. 
The (red) Lio,Lioo and Liooo are our results for a critical 
temperature parameter equal to C = 10, 100 and 1000. 



where C is the usual Riemann Zeta function. As adver- 
tised earUer, the dependence of the emissivity Q on c is 
Q (X cjj , a, large enhancement compared to the rate one 
would estimate from naive dimensional analysis. 

This expression accounts for the thermal annihilation 
of two H particles into a single neutrino pair. Since there 
are three neutrino species that all couple equally the the 
Z boson, and each of these are legitimate products of the 
reaction, the total power in neutrinos per unit volume 
emanating from a bulk of helium condensate of temper- 
ature T is 
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Implications for white dwarf cooling. The phenomeno- 
logical importance of energy loss rate through neutrino 
emission depends on the ratio of the neutrino luminos- 
ity L^ and the electromagnetic luminosity Lem- The 
comparison between the two is complicated by the fact 
that the electromagnetic emission is proportional to the 
surface area of the star while the neutrino emission is a 
bulk property. In addition, the electromagnetic radia- 
tion depends on the surface temperature while the neu- 
trino emission depends on the temperature of the inte- 
rior. Finally, the neutrino emission depends strongly on 
the density. All these factors make it necessary to use 
a model of the star structure in order to make a mean- 
ingful comparison. Since our aim here is to understand 
whether the neutrino emission mechanism we identified 
can be competitive with electromagnetic radiation as a 
cooling mechanism, we will stick to simple models of the 
star structure capable of giving trustable but not very 
accurate estimates. 

The broad features of electromagnetic radiation can 



be understood in a model where most of the star is occu- 
pied by an isothermal, highly heat conducting degenerate 
core surrounded by a thin layer of non-degenerate mate- 
rial. Assuming that radiative heat transport dominates 
over convection, using a free gas equation of state for the 
electrons and a Kramer's opacity law for the envelope 
one can derive a simple relation for the star's luminos- 
ity [21, 22]: 
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where L (L©) is the star (sun) luminosity, M {Mq) is the 
star (solar) mass and T is the nearly constant tempera- 
ture in the star interior. To be specific, we will consider 
a He WD with mass M — O.4O6M0. For such a star, the 
luminosity predicted by (15) is shown in Fig. 2 by the L^ 
solid (blue) line. A more complete model for He WDs 
including convection provides a similar relation between 
the interior temperature and the surface temperature [20], 
which then gives an estimate for the electromagnetic lu- 
minosity if we assume black body radiation. This rela- 
tion is also shown in Fig. 2, including a naive linear ex- 
trapolation (shown in dotted lines) for lower and higher 
temperatures than those computed in Ref. [20]. 

To estimate L^, we use a simple model of the star struc- 
ture given by the hydrostatic equilibrium equation: 
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where p{r) is the mass density, Gn is the Newton con- 
stant and P = P{p) is the pressure of a (relativistic) free 
degenerate Fermi gas 
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The solution of (16) and (17) provides us with the density 
profile of the star which we can fold into (14) and thence 
compute the L^y. However, (14) is valid only where the 
condensate forms. Close to the surface of the star the 
density is small and so is the condensation temperature 
Tcondi and hence one must appropriately restrict the ra- 
dial integral in L^. At the center of the star the density 
is high and the interior temperature T may be smaller 
than Tjneitj and instead of a Bose condensed phase one 
may have a crystalline phase. Consequently, to be con- 
servative we restrict the integration both at low and large 
r, so that 
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to the region where the condensate is sure to exist. The 



distances 



and 



are defined by 
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The classical crystal-melting temperature depends on the 
position through its dependence on the density: 
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but the presence of strong quantum effects in our system 
may well make the true Tmeit lower. As we will see, 
however, uncertainties in rmin are not very significant for 
Ly since the largest contributions to L^, come from large 
r. The condensation temperature also depends on the 
position through its dependence on the density. However, 
this dependence is much more uncertain. A free Bose gas 
would have a critical temperature equal to 
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with C = 1. The repulsion between the ions raises the 
critical temperature by a currently unknown amount. In 
Ref. [18] the critical temperature was found to be of the 
form in (20) with C ~ 9, but this calculation did not 
include the effects of dynamical electrons or the massless 
quasiparticle. Since the value of the constant C is the 
biggest uncertainty in our estimate we repeated our cal- 
culation for different values of C covering a wide range. 
The results of these calculation are shown in the three 
(red) lines denoted by Lio, Liqq and -/jiooo in Fig- 2, cor- 
responding respectively to the values C = 10, 100 and 
1000. For comparison we also show in Fig. 2 (the yellow 
line "Lj/") the neutrino luminosity from a fictitious 'star' 
with a constant density equal to the average density of a 
M = O.4O6M0 WD star. The slope of this last curve is 
set by the T^^ dependence in (14). The curves Lio,iioo 
and iiooo represent the competition of two effects. The 
first is the steep dependence Q ~ T^^. The second is that 
the higher the temperature is, the larger the uncondensed 
portion of the star is, and that portion does not emit neu- 
trinos. This effect is especially important because, due to 
the strong dependence of Q on the density {Q ^ p^^^^), 
most of the emission happens on the outer regions of the 
star, the first region to enter the non-condensed state as 
the temperature rises. The second factor dampens some- 
what the rapid increase expected by the r^ T^^ behavior 
until a temperature higher than Tcond{r = 0), at which 
point the condensed phase disappears from the star and 
the neutrino emission comes abruptly to an end. 



The curves in the Fig. 2 have some uncertainties that 
should be kept in mind. These are mainly due to the 
use of the T = dispersion relation for H, and to the 
reliance of our computation on the low-momentum limit 
of the dispersion relations. The latter point is an is- 
sue because away from small p, the dispersion relation 
of H is not linear, with the deviations becoming impor- 
tant for modes with energies T ~ po ^ i^A (recall from 
Ref. [12] that uia ranges between 2.1 • lO^K and 4.6- lO^K 
as ao/l goes from 30 to 50, where ag is the Bohr radius 
and I the interparticle spacing). Hence we can only re- 
ally trust our results for L where T < tua- The use of 
the T = dispersion relation for H may be problem- 
atic, since even when T <^ Tc near the center of the star, 
close to the surface one will probe regions where T ^ T^. 
thanks to the decrease in density. In general, one would 
expect that the effect of a finite temperature would be to 
lower TTiA by lowering the size of the nuclear condensate, 
which would then lower ch, increasing the luminosity. 
If the phase transition to the uncondensed phase is first 
order, as suggested by [18], then the drop in niA close 
to Tc would not be parametrically large, and such ef- 
fects should not affect our conclusions qualitatively. If 
instead the phase transition is second-order, there is an 
enhancement of the luminosity at points where v is small, 
but that enhancement should be cut off by the fact that 
the dispersion relation will be linear over a much smaller 
range of momentum. A better assessment of these finite- 
temperature effects must wait until the thermodynamics 
of this transition is better understood. 

Conclusions. The results shown in Fig. 2 indicates 
that neutrino emission from the condensed phase can be 
an important sink of energy in a He WD if and only if 
the critical temperature for condensation is quite high. 
More systematic calculations of the critical temperature 
and of the finite-T corrections to the dispersion relations 
are currently under way [23] . In the case the critical tem- 
perature turns out to be high it would be important to 
include the neutrino emission from the condensed phase 
in a realistic cooling code. Only then could one ascertain 
whether the cooling of He WD can be a smoking gun sig- 
nature for the existence of a quantum liquid phase inside 
white dwarfs. 
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